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Abstract 

There are strong theoretical indications that the minority valence quark 
distributions in the nucleon may break charge symmetry by as much as 3- 
5%. We show that a comparison of e^D deep inelastic scattering through 
the charged current, for example at HERA, could provide a direct test of this 
effect. This measurement is also be sensitive to an intrinsic component of the 
strange quark sea which leads to s 7^ s. 
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The assumption of charge symmetry is an essential element of all current, phenomeno- 
logical analyses of deep inelastic scattering data in terms of parton distributions. We recall 
that charge symmetry is the invariance of a Hamiltonian under a rotation by 180" about the 
2-axis in isospace: 



0. (1) 



It is a much more restrictive symmetry than isospin (i.e., charge independence), which 
requires [/j, Hqs] = 0, z = 1, 2, 3. (Here Jj is the generator of the z'th component of rotations 
in isospace.) A charge symmetry transformation maps n into p and d into u (and vice versa), 
so charge symmetry at the quark level requires 

rf"(x) = uP{x) ■ m"(x) = (F{x), 

dr{x)=uP{x) ; vJ'ix) =dP{x). (2) 

Charge symmetry thus reduces the number of light quark distributions to be extracted from 
data by a factor of two [0]. All parton distributions for the neutron can be expressed in 
terms of those for the proton. 

Studies of charge symmetry and charge independence in nuclear systems are very well 
developed . While the latter is quite often broken at the 5% level, the former is usually re- 
liable to 1% or better. Theoretical estimates of charge symmetry violation (CSV) in parton 
distributions were only made a few years ago. This was motivated by the successful calcu- 
lations of charge symmetric parton distributions in the nucleon In these calculations, 
parton distributions were obtained from the relation 

g(x,/x2)=M^|(X|V^+(0)|Ar)|2 5(M(l-x)-p+) . (3) 

X 

In Eq.(H), X represents a complete set of eigenstates for the residual system. The parton 
distribution qi^x^iJ?) is guaranteed to have proper support, i.e. it vanishes for x > 1. One 
can easily see that the lowest-mass intermediate states dominate the sum over states in the 
valence region and extensive investigation has shown that simple models, like the MIT bag 
||1CI| |, are capable of explaining the shapes of the charge symmetric parton distributions. 

Beginning with Sather ||^ and Londergan et al. @J^, it was realized that the dominant 
contribution to parton CSV should arise from mass differences in the eigenstates in 
Eq.(^, rather than from the quark wavefunctions. In this case, expressions for the parton 
CSV terms could be obtained from the charge symmetric parton distributions. The dominant 
effect for valence quarks turned out to be the u — d mass difference for the spectators to the 
struck quark. As a consequence, the biggest percentage effect occurred in the "minority" 
distribution d'^{x)—u'^{x), where the u—d mass difference enters twice (the valence spectators 
being (uu) and (dd), respectively). 

Suppose we define the valence CSV quantities 5dv and 5uv as: 

5dv{x) = c?y(x) — Uy^x); 

Suv{x) = u\,{x) - dv{x), (4) 
then it was found that 5dv{x) ~ —6uv{x) As the ratio dy/uy becomes small at large 



Bjorken x |T^, the relative CSV effect is much bigger in 5dy{x) / dy{x) than 5uy{x) /uy{x) 



2 



For example, Rodionov et al. found that the latter was never greater than 1%, as expected, 
whereas the former could be as large as 5-10% at intermediate values of x. Sather |^ 
suggested that 5dy{x) / dy{x) could grow to 3% - a value supported by recent work by 



Benesh and Goldman ||12[ - c.f. however, Ref. |T3 . 

In summary, it seems very likely that the minority valence distributions, dy{x) and Uy{x), 
may break charge symmetry by 3-5%. Experimental confirmation of such an unexpectedly 
large effect would be extremely important. Until now the best proposal to determine 5dy{x) 
and 5uy{x) has involved vr^D Drell-Yan 0. That analysis is somewhat complicated by 
the possibility of CSV in the quark distributions of the pion, that is a non-zero value of 
d""^ — u" . In this paper, we discuss the possibility of extracting information on CSV by 
analyzing charged current deep-inelastic scattering from an isoscalar target. In particular, 
we examine the possibilities for comparing charged-current DIS from electrons or positrons 
on deuterium, which could be performed in the future at HERA. 

At the enormous values of that can be probed at HERA, charge current (CC) weak 
interaction processes such as e~p u^X are not impossibly suppressed with respect to the 
electromagnetic process e^p —>■ e^X. In future experiments at HERA, one should have both 
e~ and e"*" beams, and current plans call for colliding beams involving heavier nuclei, such 
as D, in a few years |T^. The (e~, Ug) reaction proceeds through absorption of a W~ by the 



target partons. It couples only to the positively charged partons in the target, so that for a 
deuteron target the structure function (per nucleon) is [] 

^(x) = [u^ix) + dF{x) + + J"(x) + 2s{x) + 2c(x)] /2. (5) 

(For simplicity we denote the (e~, Ve) reaction by the charge of the virtual W absorbed by 
the target.) Eq.(|^) is true under the following conditions. We assume that the measurements 
occur at energies and momentum transfers well above heavy quark production thresholds. 
We also neglect terms in the CKM quark mixing matrix of the order of \Vuh\^ ~ \^td\^ ~ 
1 X 10~^. If we have a positron beam, then the (e"*", Ve) deep-inelastic reaction measures only 
the negatively charged partons: 

""^(x) = [rfP(x) + vF{x) + d"(x) + u"(x) + 2s(x) + 2c(x)] /2. (6) 

Taking the difference of the e"*" and e" CC cross sections one therefore has 

5F,{x) ^ F/^^^(x) - -^(x) = ^^^(^) - ^^v-(x) ^ ^^^^ _ _ ^^^^ ^ 



Eq.(0) demonstrates that any difference between these structure functions will arise 
from either parton CSV, or from differences between the s and s distributions (or c and 
c distributions) in the nucleon. Furthermore, the relation holds for all x values. Since the 
charged current weak interaction process is rare, it would be essential to have very accurate 
calibration of the electron and positron flux. Detector efficiencies would not be a major 
problem, as the signal involves prominent jets on the hadron side and large missing energy 
and momentum on the lepton side. 



We have also ignored any tiny corrections which might arise because ^ or d' ^ c^. 
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If, as is commonly assumed, s{x) = s{x) and c{x) = c(x), Eq.(^ provides a direct 
measure of the CSV in the parton distributions. On the other hand, there has been quite a 
lot of interest recently |]T5|]- |T8| in the possibility (first discussed in Ref. |T9[) that s{x)—s{x) 
might be non-zero - with at least some suggestion of experimental support for the idea ||20|| . 
We show below the estimate of s — s calculated by Melnitchouk and Malheiro [|T5|. This is 
quite dependent on the form factor at the NKA vertex which is not well known. Clearly it 
will be important to determine first whether or not d _ -pw z?^^ 6Fi) is experimentally 
non-zero. The interpretation in terms of CSV, s ^ s ot possibly both, can then be pursued 
in detail. 

In order to illustrate the size and shape of the effect expected we have evolved the values 
of 6dv and 6uv calculated in Ref. |Q - as shown at = 10 GeV^ in Fig. 1(a) of Ref. 
- to values of appropriate to HERA. The NLO QCD evolution of the distributions 
was carried out using the Mellin transformation technique. By the use of this technique one 
easily transforms the classical DGLAP equations P|,PT[] into a system of ordinary differential 
equations. This is principally done using Mellin moments. 



M^(g2) = /'dxx^-^F(x,g' 

Jo 



(8) 



to transform from x-space to complex A^-space. In order to regain the evolved distribution, 
once the evolution equations have been evaluated, the results must be transformed back to 
x-space by the inverse Mellin transform. This is performed by a contour integral in the 
complex A^-plane, 



F{x,Q' 



TT JO 



dz Im 



i4> ^—c—ze^'f p^n=c+ze^'f ^q2\^ 



(9) 



in which the contour of integration, and hence the value of c, must lie to the right of all 
singularities of in the complex A^ plane. For all the practical calculations, we have 
chosen to use the same values for z and (f) that were used in the original paper of Gliick, 
Reya and Vogt |22]. The major advantage of using the Mellin transformation method for 
computing the evolution of the moments is that it does not involve the tremendous amount 
of computing time that previous methods employed. The resulting predictions for 5uv and 
5dy at Q2 = 100,400 and 10000 GeV^ are shown in Fig. 1. 

In order to indicate the relative size of the difference in the e^D cross sections expected, 
we divide the difference F^"^^ — F^ ^ by the average of F^^^ and F^ ^ to obtain: 



R{x) = 



Fl 



W+Di 



X) 



W-Di 



F^^^ix) + Fl 



5dv{x) — 5uv{x) + 2{s{x) — s{x)) 
[u^{x) + u^{x) + d^{x) + (P'(x)] + 2(s(x) + s{x)) 

j=p,n 



R 



CSV 



[x) + Rs(x) 



(10) 



In Eq.(|TOD, Rcsv is the charge-symmetry breaking term, and Rg is the term proportional 
to s — s. In this equation we have neglected terms proportional to c and c since we expect 
them to be quite small, and in addition we know of no quantitatively reliable model for this 
quantity - see, however, Ref. for an estimate of the intrinsic charm of the proton in the 
context of the HERA anomaly. 



4 



For the combination F^^^ + ^ we evolved the CTEQ4Q parton distributions ||23 



from their starting scale (Qq = 1.6 GeV^) to the same values of used for Sdy and 6uv- 
The results for Rcsv are shown in Fig. 2, at several values of Q^. We see that the CSV 
contribution to the cross section rises from 1% at x = 0.4 to 2% at a; = 0.6. The prediction 
is not shown above x = 0.8 because the original bag model predictions were not reliable at 
very large x and the cross section is, in any case, too small in that region. 

In view of the appearance of s — s in Eq.(|^), we have taken the recent work of Melnitchouk 
and Malheiro [l^ as an indication of the possible size of the effect. These authors evaluated 
the strange quark distributions in a "meson-cloud" picture, where the s sea arises from the 
kaon cloud of the nucleon, while the s-sea arises from the spectator strange baryon, namely 
A or S. We have chosen to show their prediction |^5| with a (monopole) form factor of mass 
1 GeV, which is the largest value consistent with the latest CCFR data . Figure 3 shows 
the corresponding ratio Rg at the same values of as shown in Figs. 1 and 2. Clearly 
the order of magnitude of the effect caused by s — s is very similar to that arising from 
CSV. Nevertheless, the shapes are completely different and one would expect to be able to 
separate the two phenomena on the basis of the measured x-dependence of R. 

In conclusion, we note that although both the CSV and strangeness effects are predicted 
to be at the level of a few percent, the measurement of the difference in e^D cross sections 
at the momentum transfers typical of HERA is a very "clean" experiment. Any deviation 
from zero, at any value of x, would be extremely interesting, whether its origin lies in CSV 
or intrinsic strangeness. Since there are both electron and positron beams at HERA and 
future plans for deuteron beams, we hope that this comparison can be made in the near 
future. 
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FIG. 1. (a) The charge symmetry violation in the valence down quark distribution, (b) The 
charge symmetry violation in the valence up quark distribution. 
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FIG. 2. The charge symmetry violating ratio Rcsv defined in Eq.(lO). 
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FIG. 3. The relative contribution to the difference in e cross sections arising from a possible 
difference s — s - labelled Rg in Eq.(lO). 
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